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Abstract 

We perform a study of various anharmonic potentials using a recently developed method. We 
calculate both the wave functions and the energy eigenvalues for the ground and first excited states 
of the quartic, sextic and octic potentials with high precision, comparing the results with other 
techniques available in the literature. 
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I. INTRODUCTION 

Ordinary perturbation theory often leads to asymptotically divergent series Several 
techniques have been devised in the past in order to improve the convergence of the standard 



perturbative expansion (see, e. g., Refs. 



and references therein). 



Some of the "optimized expansions" that have been proposed are based upon the so- 
called linear delta expansion (LDE) {3,0], where the terms in the Schrodinger equation are 
rearranged in "non-perturbative" and "perturbative" pieces, with the contributions of the 
latter being minimized by introducing a suitably chosen variational parameter, that is by 
applying the principle of minimal sensitivity (PMS) 

Recently, in Ref. Q|, an improved ansatz for the choice of the parameter that allows 
for an "optimal" expansion has been proposed and tested in the case of the ground and 
first excited states of the quantum (quartic) anharmonic oscillator (AHO). Within this 
scheme calculating an observable with higher precision — that is going to higher order in 
the expansion — requires the solution of algebraic equations of increasing order. This can 
be done analytically, thus allowing one to get analytical expressions for both the energies 
and the wave functions of the AHO at any desired order. 

Given the success of the method for the quartic AHO, one would like to have a formal proof 
of the convergence of the method and to verify whether it works for more general potentials. 
In this paper we address this latter issue, extending the previous study to the sextic and 
octic AHO potentials: This is of interest for the study of realistic potentials — since they 
can often be approximated as polynomials — and also in itself, since general anharmonic 
potentials have applications in studies of nonlinear mechanics, molecular physics, quantum 



optics, nuclear physics and field theory 
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We also improve the convergence of the method for the calculation of the energy with 
respect to Ref. 10], by introducing a variational principle, that is by considering, at a given 
order, not the expression for the energy stemming from the LDE, but the matrix element 
of the hamiltonian obtained using the wave function at that order and applying the PMS 
to it: Since the resulting expression is quadratic in the wave function, one gets higher order 
contributions that turn out to substantially improve the convergence of the expansion and, 
at the same time, to satisfy a variational principle. 

The paper is organized in the following way: Section |H] contains a description of the 



method and a discussion of the computation of the energy through its expectation value and 
its relation to the variational principle. In Section ITTT1 we present our results for the quartic, 
sextic and octic potentials, comparing them with those obtained with other techniques 
available in the literature. Finally, our conclusions are presented in Section ITVl 



II. THE METHOD 

The method employed in this paper has been originally devised in Ref. and applied 
to the study of the quartic anharmonic potential. It relies on the the identification of three 
different scales in the problem, which reflect in a different behavior of the wave function: An 
asymptotic scale, which is fully determined by the potential; an intermediate scale, where 
the wave function decays exponentially, although with less strength; a short distance scale 
where the wave function is sizable. 

In this work we consider the application of our method to anharmonic potentials of 
the form V(x) = mu 2 x 2 /2 + p/(2N) x 2N and discuss the accuracy of the approximations 
obtained in this framework in the case of quartic, sextic and octic potentials (with N = 2, 3, 4 
respectively). 

The Schrodinger equation in the present case reads 



h d 2 mu 2 2 u 2N 

H + 7TTT X 



2m dx 2 2 2N 



1p n (x) = E n 1p n (x), (1) 



where p, is the anharmonic coupling, ip n {x) is the wave function of the nth excited state and 
E n its energy. 

Although one cannot find the solution of Eq. exactly, it is possible to determine the 
asymptotic behavior of ip n {x) in the region of large x (x — > oo) by substituting the ansatz 
ip n (x) cx e -7 '^ into Eq (□ )■ One obtains p = N + 1 and 7 = a/ pm/N / ((N + 1)K). 

In the spirit of Ref. [l0| we therefore write the wave function as 

Mx) = e-^\ N+l -^ 2 Ux), (2) 

where the exponential takes care of the correct behavior in the limit \x\ — > 00. Notice that 
the quadratic term in the exponential does not affect the behavior at large distances, but 
is relevant at intermediate scales. The coefficient (3 = m Q/2h is written in terms of the 



frequency Q = \/io 2 + Q 2 (Q is an arbitrary parameter introduced by hand, see below). The 



reduced wave function £ n is well-behaved and fulfills the equation [if]]: 

+ 




(3) 



We observe that no approximation has been invoked in the derivation of Eq. @. In the 
spirit of the Linear Delta Expansion (LDE) we now write Eq. J3J) as 

„N 



2mflx 




h 





m 2 Q 2 



2mE„ mQ 



N h 



<n{ x ) 



(4) 



H z " V N h 

In writing this equation we have added a parameter 6, which was not present in Eq. J5J). For 
5=1 one recovers the original equation in Eq. (j3J), while by taking 5 = one obtains the 
equation for the Hermite polynomials, corresponding to a harmonic oscillator of frequency 
Q. Notice that 5 here is used as a power counting device: As a matter of fact we will treat 
the right hand side of Eq. (@J) as a perturbation, although its size is not necessarily small, 
given the arbitrary nature of the parameter Q. An optimal choice of Q will make the right 
hand side of Eq. (jlj) small. 

We now write the following expansions: 



3=0 



nj 



E,, 



3=0 



"3 



(5) 



and substitute them in Eq. (jlj), thus generating a hierarchy of equations, corresponding to 
the different orders in 5. Such equations, which take the form of the equation for the Hermite 
polynomials in the presence of a source term, can be solved sequentially up to some finite 
order. 

Given the perturbative nature of the approach that we are using, all the results, ob- 
tained to a finite order in perturbation theory, will display a dependence upon the arbitrary 
frequency Q. However this dependence is artificial since the solution of Eq. (j3J) does not 
depend upon Q. In the framework of the LDE such dependence is minimized by applying 
the Principle of Minimal Sensitivity (PMS) jjj, i.e. by requiring that a given observable O 
(the energy, for example) be locally independent of Q: 
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FIG. 1: Comparison between the ground state (a) and first excited state (b) energy of the quartic 
oscillator calculated to order 15 directly through the expansion of Eq. (jSJ) (dashed) and the one 
calculated as {ipi\H\ifji) , with ipi at the same order (i = 0, 1), as a function of the arbitrary parameter 
VL, The dotted line is the result of Ref. Q]. We use m = 1/2, oj = 2 and \i = 8. 
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FIG. 2: Comparison between the ground state (a) and first excited state (b) energy of the quartic 
oscillator calculated through the expansion of the energy (squares) and the one calculated as 



(ipi\H\ipi) (crosses). The dotted line is the result of Ref. 



131 ] . We use m = 1/2, u> = 2 and \i 



In the present work we will enforce the PMS by using the energy of the state as the observ- 
able. We need however to make an important point: By solving the equations corresponding 
to the different perturbative orders, we obtain two different estimates of the energy of the 
solution. One corresponds to the expansion in Eq. ©, which was indeed used in [10]; the 



other corresponds to calculating the energy as the expectation value of the Hamiltonian in 
the state described by the wave function obtained in Eq. (j3J). This second choice turns out 
to provide a much better estimate of the energy: As a matter of fact, the energy calculated 
in this way not only contains contributions of higher order in 5 (the expectation value of the 
Hamiltonian is bilinear in the wave function) but it is also constrained by the variational 
principle (at least for the ground state) to lie above the exact result. Therefore the PMS 
applied to the expectation value of the Hamiltonian corresponds, for the ground state, to 
the statement of the variational principle. 

This finding is illustrated in Fig. Q where we compare the energy of the ground and first 
excited states of the quartic AHO — evaluated at order 15 — to the average value of the 
Hamiltonian — obtained using the wave function at the same order — as a function of Q. 

In Fig. on the other hand, one can see the ground and first excited state energies of 
the quartic AHO as a function of the perturbative order m, obtained after application of the 
PMS both to Ei m) of Eq. © and to ty^Htyf*). The faster convergence to the correct 
result in the latter case is quite apparent. 

It is interesting that also in the case of the excited state the approach based upon the 
minimization of the average value of the Hamiltonian seems to satisfy a variational principle. 



III. RESULTS 

Let us start by comparing the accuracy of our method with the one of Ref. which 
is based on a variational improvement of the ordinary perturbation theory that gives a 
convergent sequence of approximations. 

In Fig. |3] we display the logarithm of the percentile error for the ground state energy of 
the quartic oscillator as a function of the approximation order. The error is defined as 



100 x 



E n 



E„ 



(7) 



where E n is the "exact" value of the energy of the n-th state and E^ the approximate 
estimation to order m. The results using the present method (crosses) are compared to the 



outcome of the calculation in Ref. 14], obtained after 47 iterations, for the ground state of 
the quartic AHO (the "exact" value of Eq is taken from Ref. jl^). 




order 



FIG. 3: Logarithm of the error for the ground state energy of the quartic oscillator as a function 
of the order. The dashed line is the error of Eq. (11) of Ref. 14], obtained after 47 iterations. We 
use oj = and \i = 1. 



We now turn to compare our results with those of Ref. [13[ , where the quartic, sextic and 
octic AHO have been thoroughly analyzed with a method based on the generalized Bloch 
equation. This method calculates iteratively certain matrix elements of the wave operator 
(wave function expansion coefficients) and then uses a renormalization technique. 

In Fig. 0] we display the energy of the ground state of the quartic, sextic and octic AHO, 
as a function of the approximation order, compared with the results of Ref. ^j]. We have 
chosen the parameters of the AHO in such a way to correspond to the cases (3 = 2 (panel 
(a)) and j3 = 400 (panel (b)). Note that the values of Ref. jisj have been obtained using 
100 basis functions and several tens or hundreds of iterations. 

In Fig. |S1 we report, for the same cases of Fig. 01 the percentile error, denned as in Eq. ( J7| L 
but taking now for E n the value calculated in Ref. 13(. In other words, in Fig. El one can 
see, order by order, the discrepancy of our results from those of Ref. ^J. It is apparent that 
even for relatively low orders the agreement is quite good. 

The same considerations apply also to the excited states of the oscillators we are studying, 
as can be inferred from Figs. El and [TJ which display the same quantities (energy and error) 
of Figs. EJand El but for the first excited state. 

In Tables HI and HP we report the actual values of the energies of the ground state and of 
the first excited state calculated with our method to order 15 and compared to the results 



1.63 

[if 

1 .62 < 







1 


l 




1 




1 




1 




1 


1 






□ 












































































+ 


quartic 


























o 




























□ 








+ 


























>_ 






























o 






























+ 




























o 




























4> 


O- 


■o 


-o 


-o- 


o 


-o 


-o- 




-o 


o o 


o 






i 


+ 

i 


+ 


■ -h 

I 


■+■ 


■+■ 


--+■ 


■+■ 

1 


■ + ■ 


-t- 

1 


._! (_ 


. ._! 1 



2 4 6 8 10 12 14 16 

order 



(a) 



8 




+ 


1 








1 




i 




1 




1 








1.9 






+ 


+ 


. Hr. 


.+ . 




-+• 


■+■ 


■■+ 


■+■ 


- + • 


■+■ 


■ +■ 


■ + ■ 


_ 

■+■■-. 


7.8 








— 1— 












— i— 




— 1— 




— h- 




— 1 — 


5.22< 




o 














+ 




+ 












ttf 5-2 








..<>. 


































O- 


■<>■ 


o- 


o 




■o 


■<>■ 


■<>■ 


o 


■«■ 


■o- 


o ■ ■ ■ 


5.18 


















+ 




+ 












4.16C 


3- 


-+- 




— 1— 








^— 




—h- 




—h- 








— 1 — 


4.15 




□ 




































□ . 


-Q- 


□ ■ 


■□ 




■ o 


■ □■ 


■ D 


□ ■ 


■□ 


■D ■ 




■□■ 




4.14 


















i 




1 













I i I i I i I i I i I i I i L_ 

2 4 6 8 10 12 14 



order 

(b) 



FIG. 4: Energy of the ground state of the quartic, sextic and octic anharmonic oscillators, as 
a function of the approximation order, compared with the results of Ref. Q| (dotted). We as- 
sume m = 1/2 and u = 2; fj, = (8,12,16), respectively, in panel (a) and n = (1600,2400,3200), 
respectively, in panel (b). 
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FIG. 5: Logarithm of the error, with respect to the results of Ref. fl3j, of the energy of the ground 
state of the quartic, sextic and octic anharmonic oscillators, as a function of the approximation 
order. We assume m = 1/2 and u = 2; fj, = (8,12,16), respectively, in panel (a) and /i = 
(1600,2400,3200), respectively, in panel (b). 
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FIG. 6: Energy of the first excited state of the quartic, sextic and octic anharmonic oscillators, as 
a function of the approximation order, compared with the results of ^| (dotted). We assume m = 
1/2 and u = 2; // = (8,12,16), respectively, in panel (a) and fx = (1600,2400,3200), respectively, 
in panel (b). 
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FIG. 7: Logarithm of the error, with respect to the results of 13], of the energy of the first excited 
state of the quartic, sextic and octic anharmonic oscillators, as a function of the approximation 
order. We assume m = 1/2 and uj = 2; fi = (8,12,16), respectively, in panel (a) and /i = 
(1600,2400,3200), respectively, in panel (b). 



TABLE I: Energies of the ground state and of the first excited state calculated with our method 



to order 15 and compared to the case (3 = 2 in Table III of 



13j |. The correct digits are underlined. 





quartic 


sextic 


octic 




1.6075413 


1.6099319 


1.6413713 


E (RefJ13j) 


1.607541302 


1.609931952 


1.6413703 


Ei 


5.4757859 


5.7493494 


5.9996138 


E x (Ref.[13j) 


5.475784536 


5.749347753 


5.999607360 



TABLE II: As 


in Table [fl but for the case f3 = 


400 in Table III of [13j. 


The correct digits are 


underlined. 










quartic 


sextic 


octic 


E 


7.8618628 


5.188359 


4.1461938 


E (Ref.[13j) 


7.8618627 


5.1883589 


4.1461886 


Ei 


28.118479 


19.563145 


15.952016 


Ei (Ref.[13j) 


28.118454 


19.563130 


15.951985 



of Ref. |l3[, with AHO parameters corresponding to the cases (3 = 2 and j3 = 400 of that 
reference, respectively. 

Good accuracy is also obtained for the wave functions. In Fig. |H1 one can see the ground 
state and the first excited state wave functions of the quartic AHO as obtained in our 
approximation and by direct numerical calculation. 

The numerical calculation was performed using a Fortran program. We have tested 
the accuracy of the program by calculating the wave function of the harmonic oscil- 
lator and comparing it to the exact result. The error defined as A^(x) = 100 x 
|(-0( num )(;c) — ^( cxact )(x))/^( e:cact )(a;)| is found to be smaller than 10 -7 in the region where 
the wave function is sizable. As a result, a meaningful comparison of our approximate results 
with the numerical results is possible. 

Actually, the two curves in Fig. |H] are not distinguishable from each other at the scale of 
the figures: For this reason we display also the ratio of the approximate to the exact wave 
function. 



(a) 



(b) 



FIG. 8: Wave function for the ground state (a) and the first excited state (b) of the quartic 
anharmonic oscillator. The solid line is the numerical result, the dashed line (not visible) is the ap- 
proximation to order 15 of our method. The dotted line show the ratio R(x) = t// appr ^ (x) / ^ num - ) (x) 
of the approximate to the numerical wave function. We assume m = 1/2, to = 2 and /i = 8. 



IV. CONCLUSIONS 

In Ref. [n| a new method for the solution of the Schrodinger equation was introduced. It 
is based on the application of the linear delta expansion (optimized perturbation theory) and 
an ansatz for the wave function that explicitly takes into account its asymptotic behavior. 
The new method was applied to calculate the energies and wave functions of the ground and 
first excited state of the quartic anharmonic potential. 

In this paper we have extended the results of Ref. in two ways, first we have computed 
the energies by evaluating the expectation value of the Hamiltonian at a given order using the 
wave function obtained with the method and shown that the accuracy obtained is greater. 
Secondly, we have extended the results by computing the energies and wave functions of the 
ground and first excited states for the quartic, sextic and octic potentials. We have verified 
that the method works very well for these potentials and that in fact one is able to obtain 
high accuracy with only a few perturbative orders. We have also presented quantitative 
comparisons between our results and other methods found in the literature. 
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